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Abstract 

We study the neutral current flavour changing rare decay mode B — > Kjj within 
the framework of Standard Model, including the long distance contributions. It is 
found that these long distance contributions can be larger than the short distance 
contribution. 
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PACS:13.25.Hw, 13.40.Hq 



We have corrected a sign error in one of the terms contributing to the Irre- 
ducible (short-distance) contribution in our program. The revised estimates 
agree well with the ones given in a recent paper [12j. 



1 Introduction 



Radiative decays of the B-meson are very useful indices for testing the un- 
derlying theories of flavour changing neutral currents (FCNC) since the am- 
plitudes are very sensitive to QCD contributions as well as to possible contri- 
butions from loops with supersymmetric partner particles. Of the radiative 
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modes, the decay mode B —> X s "f has been measured experimentally and has 
also been the subject of extensive theoretical investigations. The advent of 
the B-factories will probably make possible measurement of radiative decay 
modes with two photons. Hence the basic two photon amplitude b — > S77 
has also been studied starting with the work of Lin, Liu and Yao p[j, and 
pursued further in references [2], fSj. 

The b — ► S77 amplitude falls naturally into two categories; an irreducible 
contribution and a reducible one where the second photon is attached to the 
external quark lines of the b — > S7 amplitude. The irreducible part can be 
estimated through the basic triangle graph. Evaluating the reducible part 
at the quark level then presents no problem. For the process B — > 
it is appropriate to the consider the amplitude as dominated by the quark 
level amplitude. However, when we consider an exclusive channel B —> M77 
with a specific meson M, it is more appropriate to consider the reducible 
contributions as arising out of emission of the second photon from the external 
hadron legs of the amplitude B — > rather than restricting oneself to 
quark level considerations and thereby completly neglecting the spectator 
contributions in B and M. The amplitude B — > ^77 then stands out since 
the basic one photon amplitude B —> Kj vanishes for real photons and we 
may expect the irreducible part of the amplitude to dominate. There will of 
course be the usual long distance contributions, the most important being 
the process B — > K*j followed by decay of K*, and B — ► Kt] c followed by 
two photon decay of the rj c . The corresponding 77 contribution will be small 
owing to the small 77 coupling to cc but we include it also for completeness. 
The B — > Krf amplitude however is anomalously high, probably because 
of the gluon fusion mechanism. Since rf has a sizeable branching ratio into 
2 7, the rf contribution to B — > ^77 will prove to be the largest resonance 
contribution. 

This note reports an estimate of the branching ratio for the process B — > 
i^77 taking account of the irreducible contribution as well as all significant 
resonance contributions. The earliest investigations of this process restricted 
themselves to the r] c contribution |lj ; subsequently an estimate was made of a 
possible contribution from the resonance rj' with the rj' arising through gluon 
fusion from the basic process b — * sgg and subsequently decaying into two 
photons . The rj' resonance however is very narrow and its contribution will 
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show up as a distinct narrow peak in the 77 invariant mass spectrum. This 
can experimentally be easily separated and we could restrict our analysis by 
excluding this resonance, but we include it for completeness. 

There are of course interference terms from each pair of terms. Unfortunately 
the only interference term whose sign we can determine is that between the 
irreducible term and the r] c resonance contribution. The other signs are un- 
known. We quote the largest and the smallest branching ratio, and give 
enough information for the reader to construct the other possibilities. 



2 Irreducible triangle contributions 

These arise from triangle diagrams with the photon being emitted by the 
quark loop. The effective Hamiltonian is (Ij 

n eff = -^wzawo^), (i) 

with 

0\ = {SiCj)v-A(Cjbi)v-A, 
2 = (SiCi) V - A {Cjbj)v-A, 

3 = (s i b i )v-A'E l {QjQj)v-Ai 

q 

4 = {Sibj) V - A J2(mi)v-A, 

q 

5 = (sibi) V - A J2(qjqj) v+A , (2) 

q 

6 = (sib^v-AYsiqjuW+A, 

q 

7 = Y^Sia^{m s P L + m h P R )biF^ 

and 

8 = -^SiG^ {m s P L + m h P R )T^Gl v . 

The invariant amplitude corresponding to quark level transition b — ► S77 
(with an incoming b line and an outgoing s line and the two photons being 
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emitted by the quark loop) is 
lQV2aG F V*V tb 



M 



b->s 



9tt 



u{ps) \T,A q J(m 2 q )yP L R 



(3) 



+ iB(m,K(m?,)P L + m b K(m 2 b )P R )T t 



where 



+ C(-m s L(m 2 s )P L + mbLim^P^e^apk^k^uip^^k^ih), 

Rpvp = kxyt^puxk^k^ — k2^€^p a \k^k2 + {ki.k 2 )e fll/ p (7 (k 2 — k\) , 

Tp, v = k 2fl k ll/ - (ki.k 2 )g^, 
A U = 3(C 3 -C 5 ) + (C 4 -C 6 ), 



and 



A d = ^[3(C 3 - C 5 ) + (C, - C 6 )}, 
A c = 3(d + C 3 - C 5 ) + (C 2 + C 4 - C 6 ), 
^ s = ^ = |[3(C 4 + C 3 - C 5 ) + (C 3 + C 4 - C 6 )] : 

5 = C = -i(3C 6 + C 5 ). 



(4) 



To get to .M(£> — > if 77) from the quark level amplitude (.M(& — > S77)), 
we should replace (s|r|6) by (if|r|£>) for any Dirac bilinear T. Recall that 
(s|s7 p 6|6) gives a factor u s j p Ub and similarly for sj p ^f b, S7 b, and sb, and 
also that (K\T\B) is zero for T = 7 s , 7 P 7 5 . In this way we get 
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■{K\sYb\B)Y,A q J{m 2 q )R 



[LVp 



1 

2 



+ ^{K\sb\B)\iB{m s K{ml) + m b K{m 2 ))T^ (5) 



+ C(-m s L(m 2 ) + m b L(m 2 b ))ep„ a pk^k^ 

In the above expressions we introduced the functions 

J(m 2 ) = I u (m 2 ), K(m 2 ) = 4/ u (m 2 ) - / 00 (m 2 ), L(m 2 ) = / 00 (m 2 ), 
where 



I pq (m 2 ) 



dx 



1-x 



dy 



x p y q 



m 2 — 2{k\.k 2 )xy — it 



(6) 



We use the following parametrization for the matrix elements of the quark 
vector current: 



(K\s 7 ,b\B) = Upb + Vk] 



m 



B 



m 



K 



+ 



m 



B 



m 



K 



with q = Pb — Pk = k\ + k,2- It then follows that 



(K\sb\B) = {m b -m s y l 



q,(K\s^b\B) 



,77i6 - m t 



(m 



mi)F 



BK 



Also observe that 



q p R^ P = (h + k 2 ) p R^p 



(7) 



(8) 



(9) 



For the form factors appearing above, we use the explicit numerical depen- 
dence of F(q 2 ) on q 2 given by Cheng et al JB]. 



3 Resonance contributions 



3.1 The rj c resonance 



The rj c contribution to the deacy process comes via the t-channel decay B — > 
Kt]ci with then r\ c decaying into two photons. Since this has to be added to 
other contributions, we need to be careful about the sign of the term. 

The T-matrix element for this process can be written as 



i(K 71 \T\B) = (K Vc \ - iH\B)- 2 f*Hry| " M\ri e ). (10) 

q 2 - m 2 c + imrjT^ 



Noting that in the lowest order, for any states X, Y 

(X\T\Y) = -(X\H\Y), 

we get 



( |T|B> = JKv c \T\B) {ll \T\, 
q 2 - m 2 + im Vc T^ tal 
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The amplitude {^^(\T\rj c ) is parametrized as J3j 

We can determine B Vc from the r\ c — > 77 decay rate: 



(12) 



r fe-77) = g)^-7 



(2tt) 3 2A;? (2tt) 3 2A: 2 0, 



2 7 r) 2 (5 4 (/c ?7 -/c 1 -fc 2 )|(77|r|r/ c 



(13) 



We have 



and so 



E 1(77171^ = 2\B v fml, 



spins 



r(rj c - 77) = l ^ |2m ' c - 

Next we need the B — > Kij c amplitude 

{K Vc \T\B) = -{K Vc \H eS \B) 

The relevant piece of 7{ e s for this matrix element is 
Gp 



(14) 



VcbV* 



C 2 (cb) v „ A (sc) V -A + Ci(cibj) V - A (sjCi) V -A 

Ci(cc) V ^ A (sb)v-A + C 2 (CiCj) V -A(Sjbi)v-A ■ 

The second expression is obtained by Fierz transforming the first one. Thus 

(15) 



Gj_ 
v/2 



Between \B) and \K) only the V part of (sb)v-A contributes whereas for the 
c-loop with two photons only the A part of (cc)v-A is relevant. Hence, 

n eS = ^=W2(Ci + y)(c7,75c)( S - 7 ^). (16) 

Using factorization we can write 

(Krj c \T\B) = + y)(if|57^|5)fe c |c7 M 75c|0). (17) 



Define 



q»{K\s Jfl b\B) = F (ml)(m 2 B - m\ 
We thus get 

(K Vc \T\B) = -i^V th V:j v M™l){m 2 B - m\){C x + ^) (18) 
where a dipole form of the form factor Fo(m^) is used. 

Now we fix the relative sign between B Vc and f Vc . This can be easily done 
through the anomaly equation [7j 

{ll\T\Vc) = efik^ik^Y^ (19) 
Comparing with our definition of B Vc in eq.( H2|) . we identify 

and see that the relative sign is positive. However, we emphasise that we do 
not use this theoretical result to determine the numerical value of either B n 
or f Vc , but merely use it to fix the relative sign. 

The total contribution due to r/ c resonances is thus 

M Vc = 2B v J Vc ^V tb V; s (C 1 + ^)F (ml)(ml-m 2 K ) (20) 



q 2 - ml + im Vc r v t c otal 



3.2 The T] contribution 

Analogous to r/ c , the ry-resonance will also contribute to the decay amplitude 
because of its coupling to cc channel. This contribution, .M^has exactly the 
same form as eq.(]2T)|) with the parameters B Vc , f Vc , Fo(m^), and r^ tal 
being replaced by their ^-counterparts. However in this case we cannot define 
the relative sign of the r] and any of the other components of the amplitude. 
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3.3 The rf contribution 



Like rjc and 77, the rf can also contribute via the effective Hamiltonian, eq.(15). 
However, unlike rj c and 7/, the rf has a very strong coupling to a two gluon 
state. Theoretical models for rf production in B-decay via gg-states exist 
both when the gluons arise from a basic b —> sgg process JHj and also when 
the spectator quark in B-meson emits a gluon to combine with a basic b — > sg 
process fH [S]. These estimates have considerable theoretical uncertainities. 
Fortunately, experimental data onB^ Krf exists, which is sufficient for our 
purpose since the invariant mass of the two photons will be strongly peaked 
at mi,. Calling the coupling of B l — ► K l rf (i = +,0) as F(B l K l rf), the rf 
contribution to our amplitude can be written as 

Mrf = 2B r} ,F(B i K i rf)e\» (h)e* 2 v (k 2 )e^ aP k^- — 1 (21) 

q 2 - m*, + vm^r^ 

where B v > is defined as in eq.(12) with rj c — > rf. Also, F(B l K l rf) is related to 
the decay rate T(B l — > K l rf) as: 

r(5* - K'rf) = ^—\F{B i K i rf)\ 2 X^(l i , ) (22) 
lb7rmB m B m B 

Once again the relative sign of the contribution is unknown. 



3.4 K* resonance contribution 

In this case we have the following situation: 

B\ PB ) — ►iT < + 7(Ai), i = +, 

followed by 

► J Ff i + 7 (fc 2 ), 

and the process with £4 <->• &2- 

The effective Hamiltonian for the first vertex is 

n eS = ^V cb V; s C 7 s L o [W b R F» v (23) 

and the T-matrix element T % for the process K* 1 — > ET* + 7 is 

r * = 9k*k^s(pb ~ ki)ep(k 2 )k 2a (p B ~ h^e^xs (24) 

8 



with \g l K * Kl \ to be determined from the decay rate of K* 1 — > K l +^j. Summing 
over the final spins and averaging over the initial spins we get 



Iviri 2 = 1 \ a i\ ^ m K*- m K? 

3 ' 3 2 



(25) 



so \gi\ is determined from 



y I / m ir* ~~ m K 



r = 



2 \ 3 



967T \ TJIk* 



(26) 



The matrix element of 



° 7= (£) Si<7 "^"" 



IS 



{K*\q)\0 7 \B\ P + q)) = -*(^)^[^-^(pM (27) 

e ff (g) be a ^ T g T - 2^^] 

with F* to be determined from the radiative decay B l — > K* l j. With these 
definitions the complete T-matrix element can be written as 



(if77|T|5) = .M^* 

= r^(/c 1 ,A; 2 ) + (^^ i /, /ci^/c 2 ) 



(28) 



with 



(29) 



95a' 



in 



K* 



(PB - h) 2 - m\* + miK*Tf^ al 



In the last expression we have used the antisymmetry of F^ v to write the 
terms in this particular form. 

Also note that the sign of this contribution can not be unambiguously fixed 
as was donefor the r/ c contribution. 



4 Results 



The total contribution to the process B — > ^77 is 

Mtot = M lrre d + M Vc + £k*M K * + V^r? + ^'^r?' (30) 

where t; a = ±1 are sign parameters, introduced because the signs of these 
terms are unknown. The total decay rate is given by 



dT ( 1 



77 



' 1 _s I1 + m|^ 2 4m| 



d0 sin6*|A^toi| 



2 1 2 / 2 

(31) 

where ^/s 77 is the CM. energy of the two photons while is the angle which 
the decaying B-meson makes with one of the two photons in the 77 CM. 
frame. The numerical values of various parameters used in our calculation 
are listed in Appendix. 

Figure 1 shows the spectrum of our results given as a function of the invariant 
mass of the two photons, for the case of neutral B meson decay. The figure 
shows the expected resonance peaks in the 77 spectrum at the positions of 
the 77, 7/, and 7/ c , and the K* resonance contribution is spread out across this 
projection of the Dalitz plot. We have chosen the sign parameters £ a = +1 
in all cases in the graph. Figure 2 contrasts the irreducible and the total 
contribution to the branching ratio as a function of momentum transferred 
squared over the whole kinematic range. 

The various contributions to the branching ratio are given in Table 1, with the 
maximum branching ratio (£ a = +1 for all a), and the minimum branching 
ratio (£ a = — 1 for all a). We see that we predict the branching ratio in the 
range 

1.439 x 1CT 6 < Br(B -> K77) < 1.485 x 1CT 6 (32) 

and that the largest single contribution to the branching ratio comes from 
the 7/ resonance term. In fact the resonance terms (neglecting interference) 
are far far greater than the ireducible term. 

In an attempt to enhance the relative contribution of the irreducible term we 
have calculated the partial branching ratio with a cut on ^/s 77 > m r?c + 2r ?7c ~ 
3.02 GeV. In this partial branching ratio the r) c term does dominate and is 
about four the K* contribution. We could try to further reduce the non- 
irreducible background by introducing a further cut in the Dalitz plot to stay 
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Table 1: Contributions to the B° — > fT°77 branching ratio — the total branching ratio and 
with a cut on x /sZZ 



Contribution 


Branching ratio 


With cut 




xl(T 7 


xlO- 7 


Resonance 






Vc 


3.86 


1.94 


V 


0.004 





rf 


9.69 


0.005 


K* 


1.07 


0.56 


Irreducible 


0.0195 


0.0193 


Interference 






Vc ~ I 


-0.025 


0.076 


rf — I 


=F 0.13 





Vc-K* 


=F 0.048 


=F 0.099 


rf-K* 


± 0.054 


± 0.018 


K* -I 


Tl.5 x 10" 4 


=Fl.43 x 10" 4 


Maximum BR 


14.85 


2.72 


Minimum BR 


14.39 


2.48 



above the K* peak in -Jsk^, but because of the already small branching ratio, 
we do not suggest this additional cut. 

For completeness, we also quote the analogous individual contributions to the 
branching ratio for the corresponding charged decay mode (B + — > K + jj) in 
Table 2. 

Comparing the the two branching ratios for the neutral and charged deacy 
modes, one finds that although the rf contribution is larger in the case of the 
charged mode, the overall branching fractions are almost identical, whether 
one considers the complete branching ratio, or just the part above the cut. 
This is because of the fact that the larger branching fraction of the charged 
mode B + — > K* +r y gets compensated by the smaller value of the branching 
ratio for the mode K* + — > K + j. 

From the tables above, it is clear that above the imposed cut, the r] and rf 
do not contribute significantly. But the same is not found to be true for 
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Table 2: Contributions to the B + — > K +n /~/ branching ratio — the total branching ratio and 
with a cut on ./sZZ 



Contribution 


Branching ratio 
xl(T 7 


With cut 

xl0~ 7 


Resonance 

Vc 

V 

K* 


3.85 
0.0037 
12.61 
0.9 


1.94 


0.005 
0.46 


Irreducible 


0.0196 


0.0193 


Interference 

Vc~ I 

rf — I 
Vc-K* 
rj'-K* 
K* -I 


-0.025 
=F 0.13 
=F 0.043 
± 0.067 
=Fl.6 x 10" 4 


0.076 


=F 0.082 
± 0.018 
Tl-47 x 10" 4 


Maximum BR 
Minimum BR 


17.6 
17.12 


3.61 
3.41 



rj c . In particular, the interference terms are practically unimportant, with or 
without the cut. This is essentially due to small contribution from the irre- 
ducible sector while for interference between two resonance contributions, the 
fact that the cut is practically far from the resonances results in a negligible 
contribution. 

The value obtained for Br(B — ► Kjj), is much much higher than the corre- 
sponding branching ratio for the inclusive process B — ► X s 77, where the total 
branching ratio is quoted to be (1.7 — 3.37) x 1CT 7 [3j. Even the exclusive 
branching ratio with much of the long distance resonance effects eliminated 
by the cut in J~s^ we have suggested, could exceed these values (or atleast 
be of similar size). However, it is known jH] that the inclusive rate calculation 
based on quark transition picture is not valid when M\ becomes low, which 
is true in our case. Thus we should not be too surprised by our larger result 
for just one part of the inclusive branching ratio. 

It is worth mentioning at this point that if instead of using the explicit numer- 
ical dependence of B — > K form factors, following Chen etalfTOj we assume 
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Ff K (0) = Ff K {0) = 0.341 and F^ K {q 2 ) = F ^ 0) with M po/f = 6.65 GeV 

pole 

we get an even higher value for the branching ratio. 

Observation of this decay experimentally thus is expected to provide an in- 
teresting test of the underlying theory leading to this value. In particular 
the branching ratio above a cut which removes the 77, rj' and rj c contributions 
will, ideally speaking, be a good test of the magnitude of the irreducible and 
will thus be sensitive to the effects of high mass intermediate states, even to 
physics beyond the standard model. 
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Appendix 

We give the input parameters used in the numerical calculations. 





c 2 


c 3 




c 5 


c 6 


C 7 


C 8 


-0.222 


1.09 


0.010 


-0.023 


0.007 


-0.028 


-0.301 


-0.144 



Table 3: The approximate values of C[s at /i — m& 



m b = 4.8 GeV m c = 1.5 GeV m t = 175 GeV 
m s = 0.15 GeV m u = rrid = 
m B o = 5.2792 GeV r^° tal = 4.22 x 10~ 13 GeV 
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m B + = 5.2789 GeV rg| al = 4.21 x 10~ 13 GeV 
m K o = 0.497672 GeV m K+ = 0.493677 GeV 

m Vc = 3 GeV B Vc = 2.74 x 10~ 3 GeV' 1 

r?otai = 1.3 x 10- 2 GeV f Vc = 0.35 GeV 

m v = 0.547 GeV T? otal = 1.18 x 10~ 6 GeV 

B v = 13.254 x 10~ 3 GeV' 1 f v = -2.4 x 10~ 3 GeV 

Try = 0.95778 GeV r£ tal = 0.203 x 10~ 3 GeF 

Br(B° i^V) = 4.7 x 10" 5 £r(£+ = 6.5 x 10~ 5 

m^*o = 0.896 GeV g K *o Kl = 0.384916 GeV' 1 r£* al = 50.5 MeV F l = 0.52-0.72 

m K * + = 0.8916 GeV g K *+Kj = 0.355 GeV' 1 r£ al = 50.8 MeF 
We follow the Wolfenstein parametrization of the CKM matrix with 

A = 0.8 A = 0.22 77 = 0.34 

V*& ~ 1 Vu = -AX 2 
V cb = AX 2 V cs = l-^ 
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Figure 2: The spectrum of B° — > K ^ with and without resonances. We plot the differential 
BR as a function of momentum transferred squared over the whole kinematic regime. 
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